Inflationary Scenarios with Scale-invariant Spectral Tensorial 

Index 



Cesar A. Terrero-Escalante, Eloy Ayon-Beato, Alberto A. Garcia 
Departamento de Fisica, Centro de Investigation y de Estudios Avanzados del IPN, 
Apdo. Postal 14-740, 07000, Mexico D.F., Mexico. 



Abstract 

Next-to-leading order expressions related to Stewart-Lyth inverse problem are 
used to determine the inflationary models with a tensorial power spectrum 
described by a scale-invariant spectral index. Beyond power-law inflation, 
solutions are characterized by scale-dependent scalar indices. These models 
can be used as assumption on the generation of primordial perturbations to 
test for scale dependence of scalar index at large angular scales. If such a 
dependence is detected, a nonzero contribution of gravitational waves to the 
CMB spectrum must be expected. 

PACS numbers: 98.80.Cq, 98.80.Es, 98.70.Vc 
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I. INTRODUCTION 



With measurements carried out by experiments Boomerang and Maxima- 1 |TJ observa- 
tions of cosmic microwave background (CMB) anisotropies entered a phase where different 
theories for structure formation at our Universe can be falsified. To date, inflation |||| is 
the cosmological model supported by analysis of these observations j| as favored theory. 
Almost any model that produces an accelerated expansion of the early Universe (i.e., infla- 
tion) solves the well-known set of fundamental problems faced by the Standard Cosmological 
Model, provided that the inflationary period lasts long enough (see || for detailed expla- 
nations). The simplest inflationary scenario describes the classical and quantum dynamics 
of the early Universe dominated by a single scalar field (inflaton) evolving in a nearly flat 
potential. 

Analyses such as those in Refs. |4j] and |6|-|9| consist of maximizing a likelihood function 
over the space of model parameters. A given set of parameters yields a theoretical CMB 
spectrum to be compared with observations. The precise number of parameters depends 
on the version of the analysis and commonly is as small as 6 and as large as 11. They are 
classified as inflationary or cosmological parameters depending on whether they determine 
the initial power spectrum of fluctuations or its posterior evolution. To set the initial 
conditions what it is commonly assumed is a particular kind of single scalar field model 
namely, power- law inflation ||10|| , characterized by scale- invariant scalar and tensorial spectral 
indices differing in unity each from the other. Amplitudes of tensorial perturbations are often 
neglected. The conclusion to be drawn is that, in the corresponding to observed Universe 
scales, the actual potential has a strong similarity with that of power-law inflation. No 
conclusions can be made about potential functional form in other scales, particularly at 
Planck scales. 

When the relative contribution of kinetic energy to total scalar field energy (known as 
e) may be regarded as constant then, to leading order (LO) in an expansion in terms of e, 
the scalar and tensorial indices and r, the relative contribution of primordial gravitational 
waves to CMB anisotropies, are directly related and could be taken as constants too. In 
this situation, the model that better accomplishes the job of fitting the available data is 
still power-law inflation, with some distortion taking place for the relation between spectral 
indices [[TTJ. Assumption e = const, is not realistic since the equation of state must change 
near the end of inflation in order to return to the Standard Model expansion rate. Hence, 
several authors had pointed out that some scale dependence for the scalar index may be 
expected |p~^JT3[| . Furthermore, in Ref. || it was already found currently available data to 
be compatible with scale-dependent scalar index, though this dependence seems to be highly 
constrained |T3| . If a scale-dependent scalar index is the actual case, then using a power-law 
index would affect the best-fit values of the entire set of cosmological parameters. 

In this work we would like to address the question of which is the theoretical inflationary 
potential generating the primordial fluctuations characterized by a scale-dependent scalar 
index that could give the best fit to current as well as to next generation of observations. 

Among other alternatives for deriving the inflaton potential from observations (see 
Ref. ||151 , and references therein), the Stewart-Lyth inverse problem (SLIP) was introduced 
T6| as a procedure consisting of solving a pair of non-linear differential equations derived in 



Refs. [|T7l , |l8| from next-to-leading order (NLO) algebraic Stewart-Lyth equations for spec- 
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tral indices [|19[], and determining the corresponding inflaton potential. The full power of 
this procedure can be used only when information on the functional forms of both spectral 
indices is available. Unfortunately, this information is rather difficult to be directly obtained 
from observations and, in addition, simple functional forms of the spectral indices involve 
great difficulties while solving SLIP. However, there is an alternative way of using the SLIP 
related expressions. Having independent information on the functional form of one of the 
indices, the second one can be calculated and used as prior for maximizing the likelihood 
function. The resulting maximized likelihood will approve or reject the used spectral indices 
and the corresponding inflationary potential as a reliable scenario of the Early Universe. 

Next generation of observations is expected to give information on CMB polarization 
which, in turn, is linked to r (see Ref. |HJ for a review on this subject). Hence, future 



CMB data, such as MAP and Planck [21] will provide, will determine a range of likely 
values for tensor-modes parameters, in spite of the subdominance of the gravitational- waves 
contribution to CMB anisotropies characteristic of most inflationary models. Nevertheless, 
even if a central constant value of the tensorial parameters can be estimated from these 
observations, it will be almost impossible to detect any scale- dependence for the tensorial 
spectral index. With this regard, it makes sense to state the problem as finding out which 
are the inflationary scenarios, other than power-law inflation, with scale-invariant spectral 
tensorial index. A major theoretical advantage here is that this is one of the few cases where 
the SLIP is analytically solvable. 

In this paper we propose two models that could serve to fit the inflationary perturbations, 
detecting, if there exist, scale dependence for the scalar spectral index at large angular scales 
while increasing the scale range and resolution of CMB observations. These models are 
obtained in a straightforward manner as SLIP solutions assuming a constant and almost 
negligible tensorial index. 

In next Sec. we briefly describe the theoretical frame for Stewart-Lyth calculations, the 
main features of power-law model upon which these calculations are based and current ob- 
servations are fitted, and present Stewart-Lyth algebraic NLO equations for the spectral 
indices. In Sec. Ill SLIP is rewritten using e as the basic parameter. Further, we introduce 
a criterion that allows us to determine when a given SLIP solution is consistent with assump- 
tions underlying calculations. Sec. |IV| is devoted to the main aspect of this manuscript, i.e., 
how to use NLO expressions related to SLIP to determine a theoretical potential yielding a 
spectrum of perturbations able to fit observations likely to be done in the near future. We 
summarize main results obtained in Sec.|V[ 



II. PERTURBATIONS PRODUCED BY INFLATIONARY MODELS 

Many propositions for inflaton potential can be made fulfilling the conditions for suc- 
cessful inflation (see Ref. 0] for a description of some inflationary models). Hence, the 
criterion for choosing the potential must be the agreement between theoretical predictions 
and measurements. Ultimately, testing this agreement involves the calculation of primor- 
dial perturbations spectra. These spectra are used as initial conditions for the evolution 
of perturbations which can be computed through the transfer functions (using for example 
the CMBFAST package ||), solutions being compared with current measurements of CMB 
anisotropies. The simplest scenario where that comparison can be carried out is that of 
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a single and real scalar field rolling down a potential. In this scenario, a flat Friedmann- 
Robert son- Walker universe is assumed containing a single scalar field equivalent to a perfect 
fluid with equations of motions given by 



H 2 = 1 [T + V(<j>)}, (1) 
+ 3#0 = -V'(<f>) , (2) 

where <fi is the inflaton, T = <f) 2 /2, V(<p) the inflationary potential, H = a/a the Hubble 
parameter, a the scale factor, dot and prime stand for derivatives with respect to cosmic 
time and respectively, k = Sn/m^ is the Einstein constant and mpi the Planck mass. 
In this framework, the first three slow-roll parameters were respectively defined in Ref. 



22 and can be written as 15.16 



c (^) = 3T[T + ^)r 1 = ^(^) a = ^, (3) 

^ ) = -jrr~^ = e ~vT=e = K dip> (4) 

Up to a constant, the first slow-roll parameter (0) is a measure of the relative contribution 
of kinetic energy to total field energy. By definition e > and, as it is well known and we 
shall see in details further in this manuscript, it has to be less than unity for inflation to 
proceed. This last assertion implies the potential being positive defined. We also note that 
while defining second and third slow-roll parameters it was assumed the potential to be a 
monotonically decreasing function. We shall return later to the point of assumptions behind 
the definitions and calculations to be used in this manuscript. 

Few models of inflation allow exact determination of scalar and tensorial perturbations. 



One of them is power-law [10], a scenario of inflation where: 



a{t) cx t p , H{<j>) ex exp (-J±<f\, V{<j>) cx exp (~^<P ) , (6) 



with p being a positive constant. It follows from Eqs. @, (§J), and (H) that in this case 
the slow-roll parameters are constant and equal each other, e = r] = £ = 1/p. Note that 
condition e < 1 implies p > 1. 

For perturbations produced during power-law inflation it can be shown |T{| that Ax(k) = 
10As(k)/ y/p, where As(k) and Ax(k) stand for normalized scalar and tensorial spectral 
amplitudes, and k is the wavenumber corresponding to the scale matching the Hubble radius, 
k = aH. Now, through definition of the spectral indices 

1 dlnA 2 s (k) dlnAl(k) ... 

one can see that for power-law inflation 
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n T n s {k) — 1 1 
T ~ 2 ~ 1-p 




<0. (8) 



where the inequality follows from above mentioned condition e < 1. Obviously, relation 
must be valid to any order in expansion of 1/(1 — p). Since spectral indices are constant, 
for power-law inflation, definitions ([?]) can be rewritten as 

/ k \ ns ~ 1 / k \ nT 

4W=4W r , AUk)=AUko) f , (9) 



where fc is a pivot scale usually taken to be the scale probed by COBE. From these expres- 
sions it is evident that n s — 1 and n T measure the deviation of the spectra amplitudes from 
scale- invariance. 

Up to the present, there are not general analytic expressions to calculate power spectra 



of inflationary models. Based on the power-law solution, Stewart and Lyth |19[ derived 
approximated expressions for both spectra regarding as small the deviation of higher slow- 
roll parameters from e and also deviation of e with respect to zero. These approximations 
imply the slow-roll parameters to be slowly varying in time functions. In terms of spectral 
indices, NLO expressions are 

1 - n s (Jfe) ~ 4e - 2r] + 8{C + l)e 2 - (IOC + 6)er? + 2C£ 2 , (10) 
n T {k) ~ -2e [1 + (2C + 3)e - 2(C + 1)77] , (11) 

where notation is that of Ref. |15| and C ~ —0.73 is a constant. Symbol ~ is used to 
indicate that these equations were obtained using the slow-roll expansion. Hereafter we 
shall use equal sign in our calculations, but meaning of approximation should be added 
whenever it applies. 

For a giving expression of scale factor, the Hubble parameter and potential are deter- 
mined, and then substituting definitions ©, (§), and (|) in Stewart-Lyth equations ( |l0|) 



and (|TT|), scale-dependent spectral indices are obtained. Imposing exact power-law relation 
e — V — £ — expression (§) is recovered. 

Even if it is already possible to include time dependence of e at lowest order, we shall show 
here that some interesting cases able to answer the question stated in this paper are only 
considered if NLO expressions are used. NLO expressions for the spectra have been tested 
and found to provide a high accuracy for theoretical perturbations calculations ||23|| . Even 
more, some authors have stressed that NLO expressions in terms of e will be compulsory in 
order to match analytic results with data to be obtained in the near future [|15| ,p4|j. 

For Eq. (O) constrain < remains valid, even for scale-dependent tensorial index. 
Assume n T > then, from Eq. (|TT|), we obtain inequality 1 < 1 + e < 2(C + 1) (77 — e). 
Evaluating C = —0.73 yields rj—e > 1.852, in contradiction with the approximation \rj — e\ < 
1 used to derive Eq. flTT|) . Thus, in general, 

n T (k)<0. (12) 

For the kind of models we are considering here, even if tensorial perturbations are large 
enough to have a detectable imprint in CMB anisotropics, according with definition ([/]), their 
amplitudes will vanish in wave longitudes to be probed by most gravity waves interferometers 
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||. A similar analysis for Eq. (|T0| ) yields that, although models with ns < 1 are favored, 
not particular constrain exist upon values of ns, thus allowing models with so called, blue 
spectra, i.e., ng > 1. 



III. STEWART-LYTH INVERSE PROBLEM 



In Ref. |TjJ it was shown that using definitions (|3[), (|), (|5[) and defining r = In i? 2 , 
= nr{k)/2 and A(/c) = [n s (k) — l]/2, the indices equations in terms of the first slow- 
roll parameter e and its derivatives with respect to r (e = de/dr and e = d 2 e/dr 2 ), in a 



straightforward manner, can be written as |17| , |18 



2Cee - (2C + 3)ee -e + e 2 + e + A = 0, 
2(C + l)ee - e 2 - e - 5 = . 



(13) 
(14) 



Given expressions for scale-dependent spectral indices, corresponding inflaton potential 
can be found by solving Eqs. (|13|) and (|I4"D for e, and using definitions of first slow-roll 
parameter (Q). This procedure is what we called Stewart-Lyth inverse problem [|T6 | . 

When SLIP was introduced in Ref. [16|, the inflaton potential corresponding to a given 



solution of differential equations ([13]) and (14), was determined as a parametric function of 
r. Solutions of these equations are expressed with r as an explicit function of e and are 
difficult to convert to expressions for e as explicit functions of r. Hence, it seems reasonable 
to look for a similar procedure, in terms of the first slow-roll parameter. Moreover, as we 
shall see later, using e as a parameter allows us to analyze the solution restricted to the 
interval of <j) where inflation is feasible, i.e., < e < 1. 

The expression for the potential as a function of e remains the same that in Ref. [H5 and 
is obtained from definition (|3|), 



V(e) = - (3 - e) exp [r(e)] 



(15) 



but here, instead substituting the first slow-roll parameter as function of r, we substitute 
exp[r(e)]. On the other hand, using Eq. (|3|) fllEf , 



(r) 



^2k 



dr 



+ 



(16) 



em 



where 0o is an integration constant. Changing variables, and substituting e from the first 
order equation (|14|), 



2{C+1) 



yfede 



+ 



This way, the inflaton potential is given by the parametric function, 



(17) 



V{<j>) 



V(e). 



(18) 
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Here is very important to recall that for SLIP solution (|18D to be unique, e(r) must be 
solution of both equations ([H|) and flT4f) . Need of information on the tensorial modes is also 
a conclusion stressed by Lidsey et al. [[DJ in their report about perturbative reconstruction 
of inflaton potential. With this regards, one can see that solution e = 1/p to Eqs. (|T3|) and 
(|I4|) automatically implies A — 5 — 1/(1— p) in full correspondence with relation (|8]). 



A. Consistency criterion for SLIP solutions 



Solving SLIP is not enough to state that solutions have any physical meaning. We recall 
that our calculations are based on several assumptions regarding the form of the potential 
(behavior of the inflaton as function of cosmic time) and range of slow-roll parameters. 
Hence, for SLIP to yield consistent results, conditions arising from this assumptions should 
be fulfilled by the obtained inflaton and its potential. Let us analyze these conditions in 
detail. 

To derive Eqs. fliPp and (|TTJ) and SLIP solution (0), no particular assumption was 
made about initial conditions for nor for its expected value thus, solutions of SLIP are 
not constrained to chaotic or new inflation nor to a particular energy scale. Furthermore, 
no assumption was neither made about the potential convexity so, in principle, it could 
be in any of categories related to classification given in Ref. @ and moreover, the same 
potential could have features characteristic of different categories. This is a consequence of 
dealing with NLO expressions, i.e., allowing a larger variation of slow-roll parameters during 
inflation. Particularly, SLIP solutions can be associated with a hybrid scenario of inflation 
f25fl . Next, we shall analyze possible constrains upon solutions of SLIP. Deriving Eq. ([!]) 
with respect to cosmic time and inserting Eq. (Q) it is obtained, 

T=—H. (19) 

Now, taking into account that H' = H /</>, to determine a sign for H' it is necessary to fix 
the inflaton behavior as a function of cosmic time. In this paper, it was assumed that <p > 0, 
and, correspondingly, H' < 0. Further, comparing definitions in Eq. (|3|) for the first slow-roll 
parameter and eliminating e yields, 
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T = ^H' 2 , (20) 



K 2 

and after substituting Eq. (|20|) in the equation of motion (P and deriving with respect to (f> 
we obtain, 

V = I (3 - rj) HH' , (21) 

where the definition of the second slow-roll parameter (Q) was used. 

Eqs. ([10]) and ([□]) were obtained using the slow-roll expansion, i.e, the absolute value of 
r) should be close to the value of e which, in turn, should be near zero. Therefore, (3 — rj) > 
and the sign of V is determined by that of H 1 . Hence, according to the previous assumptions, 
the potential must be a monotonically decreasing function of the inflaton. Summarizing, 
any feasible solution of SLIP should fulfill the following conditions: 
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\H0)<o' (22) 

These correspond to conditions for the inflaton to roll down the potential from lower to 
higher values of the scalar field. At least in the case we shall analyze here, it seems to be 



useful to write conditions (22) in an equivalent manner. Note that dr/dt = 2H / H < 0, then 
conditions (^) imply 

\v>o, ^ 14 >o. (23) 

Similar conditions could be derived for the case of an inflaton rolling down from higher to 
lower values by properly choosing the behavior of the inflaton as a function of cosmic time. 
For that case conditions (E3D will read, 



& > o. (24) 



ef > 

dV_ 

de 

In fact, of all the expressions used here, the sign of only affects Eq. (|ITD, the modification 
being 0(e) — > —0(e). That changes the SLIP solutions by the mirror equivalent solutions. 

On the other hand, inflation is defined as a period where scale factor grows accelerately, 
i.e, a > 0. This is equivalent to say that d(H~ l / a) / dt < 0, definition remarking that, during 
inflation, the comoving Hubble radius decrease with time. Deriving, using expression ( |i~9|) 
and definition (Q) we obtain already mentioned upper value for first slow-roll parameter, 
i.e., e < 1. By definition e > 0, then criteria (|23|) and (p4|) must be tested in the interval 
e G [0, 1). Looking at graphs of inflaton and its potential as functions of e, and plots of e(r) 
in the corresponding range, one is able to find out whether SLIP solutions will be consistent 
with underlying assumptions. Note that these conditions are sufficient for given potential 
to be inflationary but they do not ensure the inflationary epoch to be long enough. 



IV. MODELS WITH CONSTANT TENSORIAL SPECTRAL INDEX 

To settle down the cosmological parameters, the set containing the parameters that 
determine the realization of the cosmological model together with the parameters that de- 
termine the initial conditions is tuned in order to maximize a likelihood function @|| [|. 
Commonly, the initial power spectra are set to the form corresponding to power-law infla- 
tion with negligible amplitudes for primordial gravitational waves. The assumption behind 
is that during the inflationary lapse where quantum fluctuations were imprinted in scales 
currently reentering our causal Universe, slow-roll parameters behave roughly as constants. 

An important parameter is the tensor-scalar ratio of contribution to the CMB spectrum 
which can be defined as (see Ref. || for alternative definitions), 

r = 12.4^f. (25) 



In principle, the value of r can be estimated from observations of CMB polarization [|20 
To NLO, r is related with the tensorial index by |15 



S 



.4 



n T ~-2-f(l + 3e-2??). (26) 
A s 

From Eqs. fllOl) and (0), LO expressions are recovered by neglecting second order terms of 
slow-roll parameters, and from Eq. (|26|) by setting the expression within parenthesis equal 
to unity. This way, if in the desired scale range e and rj can be approximated by constant 
values, to LO, 11$, and r must be regarded as constants too. Hence, if some information 
on r is available, power-law inflation can still be the model providing the best fit to data 
given some distortion of the relation between indices. If some degree of scale dependence is 
hidden in the CMB, the error of assuming power-law will be reflected in the best-fit values 
of remaining parameters. Hence, to consider e ~ const, is a big restriction regarded as fair if 
the fit of data using a constant scalar index give a good result. Nevertheless, even now, such 
a good overall fit can also be achieved for some potentials yielding scale dependent scalar 
index, for example, running mass models || though the recent results in Ref. constrain 
this dependence to be very weak. Several authors have pointed out that observations with 
higher resolution and wider scale-range to be provided by satellites Planck and Map, and 
galaxies surveys should be able to discern a time dependent e [p3], making of power-law a 



poor assumption for the inflationary period. Therefore, question arises of which model could 
provide the best fit to upcoming data. To answer this question, an option is to find out 
if there exist models with slowly-varying tensorial and scalar spectral indices which, using 
current data, can be accurately fitted by a power-law and smoothly departs from power-law 
while making broader the range of scales or increasing the resolution of measurements. With 
this aim, one can make assumptions on the functional form for ny and look for solutions 
of SLIP and corresponding functional forms of the scalar spectral index which, in turn, 
can be compared with observations of CMB anisotropies. In general, the tensorial spectral 
index should be a scale-dependent function slowly varying close to zero thus, a well-based 
assumption for S is 

k k 

<5(ln k) = 5(\n k ) + a x In — + a 2 In 2 — H . (27) 

k k 

The first reliable observation of r to be obtained is hardly expected to detect any dependence 
on the scale Hence, from this information only a constant value for ut would be 



estimated. With this regard, in this paper, the analysis will be restricted to zero order in 
expansion fl2"7|). This approach has also the advantage of dealing with one of the few cases 
when SLIP is analytically solvable. We remark that to LO it is a nonsense to consider a 
constant tensorial index while regarding a time dependent e. Solutions presented here can 
be obtained only by using NLO expressions related to SLIP. 



A. Constant scalar and tensorial indices 

Let us start by proving that using 5 — C\ and A = C% as SLIP input (with C\ and C 2 
being some constants) just yields power-law solution. Recall that power-law inflation has a 
number of characteristic features: slow-roll parameters are constant and equal each other, 
spectral indices are constant, and spectra given by these indices are red tilted in the same 
magnitude from Harrison- Zeldovich spectra, i.e., C\ = C 2 ^ 0. To proceed with, we note 
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that in Eqs. (|T3|) and ([14]) e and its derivatives depend on r while A and 5 explicitly depend 
on k. So, we shall rewrite these equations in terms of k in the same fashion it was done in 



Ref. The conversion between the inflaton values and wavenumbers while crossing the 
Hubble radius can be done using expression (TT 



d In k k, H 



2H 



(28) 



Using Eq. ( f28"|) it is obtained, 



dink 

dr 



le- 1 
2 e 



After conversion to derivatives in term of Ink, Eqs. (13) and (14) become [26|, 



2e 



-e + 



C(e-l), 



2e 2 



e 2 - [(2C + 3)e + 1] 



1 



2e 



-e + e 2 + e + A = 0, 



(C + l)(e-l)e-e 2 -e-5 = 0, 



(29) 



(30) 
(31) 



where e = de/dlnk and e = <i 2 e/ci(ln/c) 2 . 

Differentiating Eq. (|31| ) with respect to In A; we can replace expressions for e and e ob- 
tained from this equation into Eq. ([!(]) and the following algebraic expression for e is ob- 
tained: 



e(k) 4 + Pe{kf + Q{8, 6, A)e(k) 2 + R{8, 5)e(k) + S(5) = 0, 



(32) 



where 



P = C + 2, 
Q(S, 5, A) = -{C + 1) [Ctf 

= (C + 1)C5+(2C 
S(5) = C5 2 . 



(2C + 3)5 + 2(C 
1)5, 



DA 



Roots of Eq. (|52| ) can be calculated but they are very complicated expressions not necessary 
for further analysis. They can be simply written as 



ei = €i(A,6,5), 



(33) 



where i — 1, . . . ,4. Thus, for 5 = C\ and A = Ci the solution is e = const.. Substituting 
back this solution into Eqs. (|30"D and (|3l|) we finally obtain 5 = A. Thus, as it was expected, 



power-law inflation is a trivial solution of our problem. Notice that, according with Eq. (^|) 
any inflationary potential that behaves like a decaying exponential for a range of corre- 
sponding to scales currently probed can be, in this range, approximated by a power-law. 
While this will offer good enough information on the inflationary period corresponding to 
currently observed Universe, it perhaps will hide clues about physics on higher energy scales. 
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B. Scale-dependent scalar index and constant tensorial index 



Previously it has been stressed the relevance of Eq. ( |I3D for constraining solutions of 
Eq. (|1^). As was already proved, for constant 5, a trivial SLIP solution is power-law inflation, 
corresponding to 5 = A. We shall analyze here remaining solutions involving scale-dependent 
A. The methodology used in this research to solve the problem of finding models with 
scale-invariant tensorial index but scale-dependent scalar index consists of solving Eq. fll4|) 
for e using as input a constant 5, substituting the solutions on Eq. (|13|) and deriving the 
corresponding functional forms for the inflaton potential and A. 

With the aim of obtaining an expression for the scalar spectral index as function of 
scale, let us first to express A as a function of e. It can be derived from algebraic Eq. (|3~2"D 
considering a constant S, 



A( 



2{c + iy 



e 2 + (C + 2) e + (C + 1) [(2C + 3) 5 + 1] + (2C + 1) - + C^l ■ ( 34 ) 



Now, we shall look for an expression of the comoving scale as a function of e. This relation 
is obtained by integrating Eq. (|3T|): 



l nfc(e)= (£±11 fin 



e 2 + e + 5 



de 



e 2 + e + 5 



(35) 



This way, the scalar spectral index can be expressed as a parametric function given by 

^) = { 0A l ( , e) ' 1 (36) 



The relevant inflationary parameters are the scalar and tensorial spectra amplitudes 
evolved through the transfer functions into the CMB anisotropies spectrum. Having Eq. 
a parametric expression for the scalar amplitudes As(k) can be readily derived using defini- 
tion (0) and Eq.flST]): 



A s (k) 




(37) 



with As(e) given by 



A s {e) = A exp 



(C+l) /A(e) 



e + 6 



-de 



(38) 



where A is an integration constant, A(e) is given by expression (|34]) and the scale as function 
of e by Eq. @. 

As it was discussed in Sec. |I| while analyzing constrain given by Eq. ([12]), we have that 
possible values of the tensorial index are restricted to the interval S < 0. With this regard, 
in Sees. |I V B 1| and [IV B 2| we shall solve the problem of finding models with scale-invariant 
tensorial index and scale-dependent scalar index for S = and 6 < respectively. 
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FIG. 1. Solution of the first order equation for 5 = 0. 

1. Null tensorial index, 5 = 

We begin by analyzing the case 5 = which, being easy to explicitly solve, is a good 
example of using SLIP to test the reliability of given functional form for tensorial index and 
of corresponding analysis of solutions. However, we shall see that the results fail to match 
current observations. 

Explicit integration of Eq. ([14]) with 6 = yields, 



exp(r - r ) = (e + 1 



.2(0+1) 



(39) 



where r is the integration constant. This solution, first reported in Ref. |TB[| , is plotted in 
Fig. |l|. Corresponding to 5 = expression for <j)(e) is obtained by straightforward integration 
ofEq. (fTJ, 



4(C + 1) 
V2k 



arctan 



(40) 



Expression for V(e) is obtained by substituting solution (|39|) in Eq. (|I~5|). 

V(e) = V (3-e)(e + l) 2 ^, 

where V = exp(r )//€. Separating e in the expression for and substituting in Eq. 
finally yields for the inflaton potential, 



(41) 



v{<p) = v - 



3 - tan 2 






cos 4 ( c+1 ) 







(42) 



Plot of this solution for range corresponding to < e < 1 is presented in Fig. 0. 

Looking at expression (|42|) and Fig. one can conclude that exist an interval of <fi where 
the potential is not consistent with the assumption <ft > 0, i.e., where the potential increases 
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FIG. 2 


. SLIP solution for 5 = 0. 





with the inflaton value. One can wonder whether there is another sector of this potential 
where inflation can take place. Indeed, such a sector can exist but the primordial fluctuations 
generated during the corresponding inflaton rolling down could be different to those assumed 
as SLIP input. Note that the range of used to plot solution ( |3SD was possible to determine 
using Eq. ( f40|) which describe the inflaton as a function of the first slow-roll parameter. The 
valid range of values for e determine the corresponding inflaton values. 

Let us look at plots of <fi and V as functions of e, Fig. 0. Regarding criterion (|23|) 





FIG. 3. Inflaton and its potential as functions of the first slow-roll parameter for 5 = 0. 



and Figs. [I], and |3], it can be concluded that solution is consistent only for e £ [0, ej$] with 
6q = (6C + 5)/(2C + 3) ~ 0.4. This way, the correct functional form for the potential with 
5 = is that in Fig. |j. Observe that for this case the potential curvature undergoes changes. 
Such a behavior for the potential is obtained as a SLIP solution due to the use of NLO 
expressions. 
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After integrating Eq. (|31|), the comoving number fcasa function of e is given by 



(43) 



with fc being the integration constant. From here, e(k) is obtained, 

1/(C+1) 



e(*) 




1 



fc . 



V(c+i) 



k_ 

k Q 



V(c+i) 



(44) 



where the root was chosen corresponding to < e < 1. 

Substituting expression fl44[) in Eq. (p4|) with 5 = 0, the corresponding expression for 
A(fc) is, 



A(fc) 




i/(c+i) 




fco. 



1/(C+1) 



x < 



\ V(c+i) 



+ 2C 



'A' 



V(c+i) 



\ V(c+i) 



(45) 



Plot for A(/c) and 5 = is presented in the left part of Fig. [5]. The range of scales plotted 
was chosen to allow observation of features at lowest scales. The logarithmic plot of a proper 
range of scales is presented in the right part of this figure, where in the vertical axis is plotted 
the logarithm of squared scalar amplitudes given by, 
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V(c+i) 



^xl/(C+l) 



-4 



14 



2- 

A 












1.96- 












1.92- 












1.88- 

1 












2(1+C) 




1.84 





10 


20 


30 


40 k/k 50 




FIG. 5. A as a function of the comoving number k and ln^4| as a function of the lnfc for 5 = 0. 



+ 



8 (C + 1) 




i/(C+i) 





V(c+i) 



2 + 4C 



\ 





i/(c+i) 



(46) 



The value for A must be chosen taking into account the observational constrain A 2 S ~ 10 -5 
given by COBE measurements. 

From these figures it is observable that the scalar index could be regarded as scale 
independent for relevant to measurements scales. The corresponding constant value is ris = 
(C + 2)/(C + 1) ~ 4.7, which is too far from values allowed by theory and experiments. 
It seems to be not possible an inflationary model to exist such that the tensorial spectrum 
generated in its framework will be of Harrison-Zeldovich type and the corresponding scalar 
index will be scale-dependent. Therefore, to this order and considering Eqs. ( p6"D and (|TTD 
with 5 = 0, it can be concluded that any model with a scale-dependent scalar index matching 
observations will give a nonzero tensorial contribution to the CMB spectrum. 

Note that, while applying criterion ( |2"4D to the mirror image with respect to the ordinate 
axis of solution (f42D, the same result is obtained. 



2. SLIP solution for negative tensorial index, 5 < 

We shall find the inflationary potentials producing perturbations that, to next-to-leading 
order, are characterized by constant and negative tensorial index and scale-dependent scalar 
index. Since for 5 < 0, SLIP equations are not explicitly solvable in terms of 0, we must 
look for a parametric expression for the inflaton potential. 

The solution of Eq. (|TJ) for 5 < is |I§ 



15 



exp(r - T ) 



e 2 + e + 5 



c+i 



2e + 1 + v / T^45 



C+l 



2e + 1 - VI -45 

The three branches corresponding to this expression are plotted in Fig. 



(47) 

It is observed 




FIG. 6. Solution of the first order equation for 5 < 0. Each curve branch that goes from left to 
right is a solution with different initial values. 



that, along with stationary solution 



2 2 



(48) 



solutions can increase unbounded or decrease bounded by value e = 0. Regarding that r 
and cosmic time have opposite signs, there is some time interval where e ~ e+, i.e., the 
solution asymptotically behaves like that of power-law. After integrating Eq. (|TT| ) for 5 < 
and inserting solution ( f47D in expression (|15|), the parametric potential is given by 



v{<p) 



2(C+1) 1 



1 + VI - 4<5arct 



an 



/2c 



+ | J-l + VT^^ln 



+ 



(49) 



t V(e) = V (3 - e) |e 2 + e + 6 



c+i 



2e+l+ v / T^4cT 



2e+l- v / T^4cT 



C+l 



The inflaton and corresponding potential as functions of e are respectively plotted in Fig. [7]. 
Note that the graph of the potential has also a maximum around 



(6C + 5) + J{6C + 5) 2 - 45{2C + 3) 



(4C + 6) 



0.4 



for S = —0.01, not shown in the figure in order to observe the details for small values of e. 
Hence, the analysis for 5 < should be done for three intervals of e, namely, I\ = [0,e + ), 



; +; 1 



and 
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FIG. 7. Inflaton and its potential as functions of the first slow-roll parameter for 5 = —0.01. 



Making use of criterion (|23|), consistent SLIP solutions are determined to exist in the 
intervals 1% and 12- Corresponding plots are presented in Fig. Here, the value of the 
constant V could be chosen to make the potential flat enough for conditions of successful 
inflation to be satisfied. For e G -Z3, the solution fails to fulfill criterion (|23|). 

The main difference between both solutions is the curvature of the inflaton potential near 
the origin, but is not a trivial one. For the same parametric expression of the potential we 
have two rather different realization of inflation depending not only on initial conditions for 
4> but also for <fi (i.e., depending on the initial conditions for e). In the case corresponding 
to the left graph on Fig. || the term in Eq. (0) given by the first derivative of the potential 
will generically dominate during the inflationary epoch meanwhile, in the remaining case, 
the evolution of <fi will be dominated by the friction term due to Universe expansion, 3H<p. 

For 5 < 0, after integration of Eq. (|35| ) we have that 



k = kn 



e 2 + e + 5 



3(P+1) 



C+l 

2 


2e+l + Vl 


-45 


2^1-44 




2e + 1 - y/T 


-45 





(50) 



where k$ is the integration constant. With A(e) given by Eq. (|34|) , parametric plots for 
A(k) corresponding to e G I\ and e G I2 are presented in Fig. |9|. Again, the scale range was 
chosen to allow details observation at lowest scales. Observe that in both cases the present 
error for ns (e.g., n$ = 0.99 ± 0.09 in Refs. Q]) can mask the scale dependence at large k^ 1 . 
After integrating Eq. (|38|) for 5 < 0, the logarithm of squared scalar amplitudes is given 

by 



In 



5C-1-C 

c^Ti 

35 (C + l) 
VI -45 



hie + 5 (C + l) In 



e 2 + e + 5 



+ 



In 



2e + 1 + VI -4:5 



2e+ 1 - VI -45 



+ 



e 3 + 2Ce 2 + 2C5 
2{C + l)e 



(51) 
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FIG. 8. Consistent SLIP solutions for e G I\ and e G I2, respectively, with 5 = —0.01. 



As observed in Fig [10], where the parametric plots of scalar amplitudes for e E I\ and eG/ 2 
is presented in the same figure, differences are almost impossible to note when the full range 
of scales is considered. Differences arise at large angular scales which could be out of reach 
for measurements. It means that, from the observational point of view, there could be not 
differences between these two realizations of inflation if the scales where differences arise are 
not probed or the resolution is not high enough to detect the scale dependence. These scales 
are precisely those where higher energies physics could leave an imprint. 

Recalling the behavior of solution (^) for small r, i.e., large t (see Fig. ||) and taking 
the limit e — > e + of Eqs. ( p£f ) and ( [BHD it is found out that for small fc -1 , A ~ 5 (see also 
Fig. |9|). That the scale interval where this behavior is observed corresponds to sufficient 
large number of e-folds to solve the Standard Model problems is provided by the asymptotic 
behavior of e near the value e+ . 

These models have the desired feature of an almost negligible 5 = const, and the scalar 
index being nearly constant in a wide range of scales, with the possibility of choosing values 
that can accurately match current observations. In fact, the value 5 = —0.01 chosen in the 
figures of this section corresponds (in the limit e — > e + ) to 11$ = 0.98 compatible with values 
given in Refs. ||. Any other constant value for ns arising from analysis similar to that 
of Refs. H,[6|-3, except for blue spectra, can also be fitted with models given by Eq. ([191). 
The weak scale-dependence of n$ obtained is in good agreement with results in Ref. |14 . 
Furthermore, the above mentioned value of 5 approximately corresponds to r ~ 0.12 which 
is greater than r* = 0.1, the value given in Ref. as the lower limit for r to be detectable 
with a 95-percent confidence regarding the error reported in Ref. M as the estimate for 
Planck measurements. Lower values for r can be appropriately taken into account. 

Because models given by Eq. (f49|) do not have a graceful exit to the Standard Model 
stage of Universe evolution (e converges to e + not to 1), 4> must be regarded here as the 
dominant scalar field in a hybrid scenario with e + being the value corresponding to the 
critical value of near which the false vacuum becomes unstable and the multiple scalar 
fields roll to the true potential minimum ||. 
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FIG. 9. A as a function of the comoving number for e € I\ and e 6 /2> respectively, and 
<5 = -0.01. 



Finally, all of the statements done in this section with regards to solution fl49|) are also 
valid after changing 0(e) by —0(e) and using criterion (p4j ). 



V. CONCLUSIONS 

We presented a version of Stewart-Lyth inverse problem using the first slow-roll param- 
eter as the basic variable in the procedure of finding the inflaton potential. That allows 
us to analyze the solutions in the range of this parameter where inflation is feasible. A 
criterion was introduced to check for solutions consistent with the assumptions underlying 
the derivation of the Stewart-Lyth inverse problem equations. 

It was shown that expressions related to Stewart-Lyth inverse problem can be used to 
determine inflationary models corresponding to given observations. We proved that power- 
law inflation is a trivial solution of this problem when constant spectral indices are used as 
input in related equations. Next-to-leading-order in the slow-roll expansion makes possible 
to consider more general scenarios where slow-roll parameters can slowly vary with time. In 
a near future, these scenarios could be more realistic than common assumption regarding 
slow-roll parameters as constants during inflation. 

Looking for a potential generating the primordial perturbations able to grow into CMB 
anisotropies and matching current and future observations, we solved the Stewart-Lyth 
inverse problem with constant tensorial index as input. Inflationary models were found 
which, unlike power-law inflation, yield scalar modes characterized by a scale-dependent 
index. For negative tensorial index, solutions were given as an expression depending on the 
first slow-roll parameter. 

The special case of a Harrison- Zeldovich spectrum of tensorial perturbations, i.e., con- 
stant amplitudes of gravitational waves, is ruled out by comparison of our results with 
current observations. It means that it seems to not exist inflationary models with scale de- 
pendent scalar index and null tensorial index. Hence, for any model exhibiting some degree 
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(upper branch). 

of scale dependence of the index of curvature perturbations it must be expected a nonzero 
contribution of primordial gravitational waves to the amplitudes of the CMB spectrum. 

Potentials obtained for strictly negative tensorial index fulfill the conditions for successful 
inflation. Evolution of the scalar field given by these potentials can be considered as the 
dynamical element in a hybrid inflation scenario, the value of the inflaton corresponding to 
power-law solution acting like the instability value for the false vacuum. 

These models can be used as assumption on the origin of primordial perturbations to 
test for scale dependence of the scalar index. Using them to fit inflationary perturbations 
is almost as easy as using power-law inflation. Only one parameter, the constant tensorial 
index, need to be fitted. If CMB polarization fails to give a value of the tensor-scalar ratio 
greater than the threshold value r* = 0.1, then the tensorial index must be fitted in the 
interval (—0.02, 0). Otherwise, if some value for the tensor-scalar ratio is measured, then an 
approximated value for the tensorial index can be estimated to serve like pivot value for the 
fitting procedure. 

We would like to stress that if any of the potentials here presented makes possible to 
reach an overall good fit for CMB anisotropies and to detect scale dependence of the scalar 
index from the next generation of observations, the conclusion to be drawn is that the actual 
inflaton potential is similar in the probed scale to the used one. In turn, if the quality of the 
overall fit is not improved compared to the result obtained using power-law inflation, either 
the scale dependence of the scalar index is practically negligible, or the information on the 
scale dependence of the tensorial index is fundamental in order to account for the features 
in the CMB anisotropies. 

The spectra of scalar perturbations produced by these models differ from those of power- 
law inflation at scales corresponding to earlier times in the Universe evolution. Thus, if while 
increasing the quality of observations a good fit using any of these potentials is achieved, it 
could give some hints about physics taking place at very high energies. 
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